Early radiative properties of the developments of time symmetric, conformally flat 

initial data 
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Using a representation of spatial infinity based in the properties of conformal geodesies, the first 
terms of an expansion for the Bondi mass for the development of time symmetric, conformally flat 
initial data are calculated. As it is to be expected, the Bondi mass agrees with the ADM at the sets 
where null infinity "touches" spatial infinity. The second term in the expansion is proportional to 
the sum of the squared norms of the Newman-Penrose constants of the spacetime. In base of this 
result it is argued that these constants may provide a measure of the incoming radiation contained 
in the spacetime. This is illustrated by means of the Misner and Brill-Lindquist data sets. 
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Introduction. The study of the radiative proper- 
ties of dynamical spacetimes describing isolated systems 
adquires these days more relevance in view of the begin- 
ning of operations of several interferometric gravitational 
wave detectors (LIGO, GEO 600). One of the crucial ra- 
diative properties to be analysed is the mass loss due 
to the radiative processes. Strictly speaking, the radia- 
tive properties are global features of the spacetimes, and 
thus they can only be unambiguosly defined at infinity. 
Penrose [0 has put forward the idea of the description 
of spacetimes describing isolated system by means of an 
unphysical spacetime (M,g^ v ) obtained from the origi- 
nal (physical) spacetime (M,g M „) through a conformal 
rescaling: = Q 2 ])^. It is within this framework that 
notions like the mass loss due to radiative processes can 
be rigorously defined. The mass loss notion is captured in 
the form of the Bondi mass which is given as the integral 
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over cuts of null infinity. The functions ^2 an d &o are 
the leading terms in of the component ip2 (C(0 3 ))of the 
Weyl tensor and the Newman-Penrose spin coefficient a 
(C(0 2 )) respectively in a so-called Bondi system. 

To the author's knowledge, the formula (|l|) has never 
been evaluated on a realistic spacetime. Furthermore, 
the dependence of the Bondi mass on the initial data (as- 
suming that the spacetime arises as the development of 
some suitable Cauchy data) whose development gives rise 
to (weakly) asymptotically simple spacetimes has never 
been analysed. It is the purpose of this letter to show 
that by means of the conformal field equations and an 
adequate gauge choice it is possible to obtain expansions 
of the gravitational field in which the dependence on ini- 
tial data can be read directly. In particular, this ap- 
proach enables us to calculate an expansion for the Bondi 
mass which shows a remarkable, yet misterious, connec- 
tion with another global property of the spacetime: the 
Newman-Penrose constants — see equation (O). Our 



FIG. 1: The standard representation of spacetime in the 
neighbourhood of spatial infinity. Spacetime corresponds to 
the region outside the cone. 



analysis will be restricted to time symmetric, conformally 
flat initial data. 

Spacetime close to spatial infinity. The natural arena 
for a discussion of the early radiative properties of space- 
time is the neighbourhood of i and <y + . In particular, it 
has been shown that under by fiat assumptions that the 
Bondi mass tends to the ADM mass as one approaches 
spatial infinity along future null infinity |J. In order to 
discuss the properties of spacetime in this region, one 
would like to make use of the Friedrich conformal field 
equations. These provide a regular system of partial dif- 
ferential equations for the field quantities in the unphys- 
ical spacetime. However, in the way they were first given 
(see for example 0) they are not too well suited for an 
analysis near spatial infinity. It can be shown that for 
the Weyl tensor one has ipabcd = (^(m/p 3 ), where p is a 
geodesic distance to spatial infinity on an initial Cauchy 
hypersurface S, and m is the ADM mass. Thus, the 
standard representation of spatial infinity as a point (see 
figure [j]) is too narrow to discuss fully the consequences of 
the conformal field equations. Using a extended version 
H of the conformal field equations which allows to con- 
sider Weyl connections and a gauge based on the prop- 
erties of conformal geodesies, it is possible to arrive to 
a new representation of spatial infinity which allows to 
formulate an initial value problem in the neighbourhood 
of spatial infinity such that: (i) the equations and ini- 
tial data arc regular, (ii) the location of null infinity is 
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FIG. 2: The cylinder at spatial infinity. Again, spacetime 
corresponds to the region outside the cones and the cylinder. 

known a priori via a conformal factor O that can be read- 
off from the initial data once the constraint equations 
have been solved ||. The unknowns in this setting are 
the frame (c p b , p, = 0, . . . , 3), the connection (T a b c d), the 
Ricci curvature (Q a bcd) and the Weyl curvature {(j) a bcd) of 
the unphysical spacetime, a,b,c,d = 0, 1 being spinorial 
indices. In this conformal geodesic gauge representation 
of spatial infinity, the point i is blown up to obtain a 
cylinder I — S 2 x (—1, 1) = {p = 0, |r| < 1} (see figure 
||). The sets = {p = 0, r = ±1} where null infinity 
touches spatial infinity are of a special nature and will 
be further discussed later. 

In the gauge given by the conformal geodesies, the 
equations for the unknowns v — (c^ fc , T abc( i, Qabcd) are 
of the form: 

d T v = Kv + Q(v,v)+L(/), (2) 

where K, Q denote respectively linear and quadratic 
functions with constant coefficients, and L denotes a lin- 
ear function with coefficients depending on the coordi- 
nates and such that L\j = 0. Differentiating formally 
with respect to p p times and evaluating on / one ob- 
tains, 

p-i 

+ (Q(v® ,v<*-fi) +£°V P ~ J) ) + L (p) (o) ,(3) 
j=i 

where denotes the operation of taking the p-th deriva- 
tive with respect to p and evaluating at p = 0, and 
<t> = (0o, 0i, 02, 03, <M- £ (p) 7^ for p > 1. The equa- 
tions for the components of the Weyl tensor are slightly 
more complicated, 

A"0„0 = £?(r)0, (4) 

where A^, p — 0, . . . ,3 are 5x5 matrices depending on 
c ab such that, 

A !/ = \/2diag(l + r, 1, 1,1,1- r), A x |/ = 0, (5) 

and B(T) is a matrix value function of the connection co- 
efficients. From the second equality one sees that the p 



derivatives disappear from the field equations upon eval- 
uation on /. Therefore, we call / a total characteristic of 
the system (||)-(^). Differentiating formally with respect 
to p p times and evaluating on I one obtains, 

A°cV0 (p) + A c d c ^ [p) = BiTfld)^ 

+ E ) { B ^abcM (P - 3) - (^)«^0^} ,(6) 

where C = 2,3. The system composed by (||) and @ for 
a fixed p are called the transport equations of order p on 
the cylinder at spatial infinity. The transport equations 
allow us to obtain expansions for the field equations of 
the form, 

■* — ' p\ 

where u&> = (r, x c ). In the case of time symmetric 
initial data, the solutions of the system (|^) with p = 
coincide with the p = solutions of Minkowski space- 
time. For p > 1, the equations are linear. The equations 
(||) are essentially ode's which allow us to calculate 
if itW for < k < p are known. Once this has been 
done, one can calculate 0( p '. Note however, that because 
of the first equality in (||), the system (||) is singular 
at the sets I ± = {p = } t = ±1} where null infinity 
"touches" spatial infinity. This phenomenon can be un- 
derstood by recalling that while is a characteristic of 
the conformal field equations, / is a total characteristic. 
Thus, a degeneracy occurs at 1^. 

In |(| a first analysis of the properties of the trans- 
port equations under the assumption of time symmetry 
of the initial data (i.e. vanishing extrinsic curvature of 
the initial hypersurface) was carried out. It is found that 
a necessary condition for the transport equations to have 
smooth solutions at the sets I ± is, 

D(a sbs ...D aibl b abcd) {i°) = 0, s = 0,l, ... , (8) 

where D ab denotes covariant derivatives with respect to 
the initial 3-metric, and b a bcd is the Cotton-Bach spinor. 
The vanishing of the Cotton-Bach spinor characterises 
locally conformally flat initial data. If the regularity con- 
dition (||) does not hold for a certain p it can be seen that 
the transport equations develop logarithmic singularities 
at 7 ± . It is conjectured that the regularity condition is 
also a sufficient condition for the solutions of the trans- 
port equations with time symmetric initial data to be 
smooth at I ± . This however remains to be proved. The 
generalisation of the condition (||) for non-time symmet- 
ric is also unknown. The importance of the transport 
equations (|3|) and (^) lies in the fact that they allow us 
to relate directly quantities defined on — the Bondi 
mass in particular — with the initial data from which the 
asymptotically simple spacetime arises. 
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Conformally flat initial data. The regularity condi- 
tion (|J) is satisfied in a trivial way for time symmetric, 
conformally flat initial data. This fact allow us to con- 
jecture that the developments of this class of initial data 
are candidates for spacctimes with a smooth null infin- 
ity. Conformally flat initial data are of great interest as 
all the time symmetric initial data sets for head-on black 
hole collisions currently in use in numerical simulations 
are of this type. We will restrict the following discussion 
to them. Being time symmetric, the initial 3-geometries 
we want to consider are completely described in terms of 
the following line element, 
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(9) 



with x satisfying the Laplace equation. It is customary 
to write, x = 1 + pW, where p is the geodesic distance 
from i. From the Yamabe equation one has that, 
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where m is the total (ADM) mass, and w p ^p,n-k 
are complex constants satisfying the reality condition 
w Pt 2 P ,2p-k = w Pt 2p,k- The functions Y n>m are the stan- 
dard spherical harmonics. Although the system of trans- 
port equations is of large size (50 x 50), its properties 
make it quite amenable to a treatment by means of a 
computer algebra system. A series of scripts in the com- 
puter algebra system Maple V have been written in order 
to solve the transport equations (|J) and (||) up to p = 5 
for time symmetric conformally flat initial data. The so- 
lutions of the transport equations up to this order are 
polynomial. This is in itself an outstanding feature, and 
it is believed to occur at all higher orders. Yet a proof 
of this fact lies far in the future. This adds further sup- 
port to the conjecture that time symmetric, conformally 
flat initial data sets yield developments with a smooth 
null infinity. The details of the implementation of the 
transport equation solver and the detailed structure of 
the solutions for generic time symmetric initial data will 
be presented elsewhere. It should be pointed out that the 
transport equations are extremely sensitive to any small 
change in their structure. Any change would inmediately 
give rise to non-smooth solutions containing logarithms. 
The absense of such kind of solutions in the expansions 
up to order p — 5 is in itself a guarantee of the correctness 
of our results. 

Calculating the Bondi mass near spatial infinity. A 
natural application of the expansions is to calculate the 
first orders of the Bondi mass in aregion of J? close to 
spatial infinity. The formula ([!]) cannot be used directly 
as it stands for it is given in the so-called Bondi gauge, 
while the gauge used in the representation of spatial in- 
finity is based in the properties of conformal geodesies. 



The connection between the two gauges has been given 
in . Knowledge of the solutions to the transport equa- 
tions up to p — 5 order permits the calculation of the first 
2 terms of the expansion for the Bondi mass. A lengthy 
calculation yields, 
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where u is a standard Bondi retarded time (it — > 
one approaches to / + ) and, 

G-l = (»™«2;4,3 - 2\/l2u;i ; 2,2Wl;2a) > 



(11) 
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(mu7 2; 4,i - 2\/l2u;i ; 2,oWi ;2 ,i) , 



are, up to an irrelevant numerical factor, the Newman- 
Penrose (NP) constants |8| of the development expressed 
in terms of the initial data as given by Friedrich and 
Kannar. Note that close to / + the second term in ( |Tl| ) 
is negative, in agreement with the fact that the Bondi 
mass is non-increasing. The formula (11) is rather unex- 
pected. Before having performed the calculation, there 
was principle, no reason for this to be the case. This ev- 
idenciates the power conveyed in our approach, and the 
kind of information that it is possible to extract through 
it. The NP constants are defined as integrals over cuts of 
of the derivative with respect of a parameter of the 
generators of light cones escaping to null infinity of the 
-00 component of the Weyl tensor. They are absolutely 
conserved in the sense that they have the same value for 
any cut of J^ + . For spacetimes admitting a conformal 
compactification that includes the point i + (future time 
infinity), the NP constants are equal to the value of the 



Weyl tensor precisely at 



Qrt. Thus, the NP con- 



stants contain global information of the spacetime that 
can be read directly from the initial data! The NP con- 
stants have also been found playing a role in determining 
the decay rate of self- gravitating waves in a Schwarzschild 
background ]Tl|: waves with vanishing NP constants are 
found to decay faster. On the same lines, an analysis of 
some particular boost-rotation symmetric solutions near 
i + has shown, remarkably, that the squared norm of the 
NP constants appear again in the leading term of the 
time derivative of the Bondi mass (late radiated power) 
p"2| . There has been some discussion concerning the idea 
of the NP giving a measure of the spurious radiation con- 
tained in the spacetime in general, and via formula ( p"l| ) 
on the initial Cauchy hypersurface in particular. The 
idea being that the first radiation reaching J? + as the 
system begins to evolve can not be truthly attributed to 
the processes occurring in the interior of the spacetime (a 
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black hole collision, for example) as it requires some time 
to reach the asymptotic region. Formula (|lTJ) is local, 
and thus valid for an arbitrarily small neighbourhood of 
spatial infinity. It will register radiation that is already 
in the asymptotic region and that had to be put in in 
order to construct the initial data. Note that formula 
( pi] ) implies that the NP constants of static spacetimes 
(the Schwarzschild solution in particular) arising from 
conformally flat data must vanish. This latter fact can 
also be independently verified using the Friedrich-Kannar 
formula for the NP constants on static initial data. Sim- 
ilarly, the recent examples of non-trivial asymptotically 
simple spacetimes given by Chrusciel and Delay |fl3|| have 
also vanishing NP constants, and accordingly a very spe- 
cial radiation content. 

The formula ([ll]) has only been calculated for devel- 
opments of time symmetric, conformally flat initial data 
sets. It is expected that a similar result would hold (mod- 
ulo some more involved calculations) for general time 
symmetric data. For initial data with non-vanishing sec- 
ond fundamental form the situation is less clear. The 
study expansions of the gravitational field obtained from 
the transport equations (0) and (||) were first carried with 
the purpose of laying the first steps towards a global exis- 
tence proof of non-trivial spacetimes with a smooth struc- 
ture at null infinity in general, and the verification of the 
conjecture regarding the regularity condition (0) in par- 
ticular. However, as (11) evidenciates they also happen 



■ Brill-Lindquist 
• Misner data 



to be a powerful tool in order to discuss the physics con- 
tained in these spacetimes. 

Comparing the Misner and Brill-Lindquist initial data 
sets. As an application of our result, one can perform 
a comparison of the NP constants of the development 
of physically equivalent Brill-Lindquist (BL) and Misner 
initial data. The interest of this example lies in the fact 
that both initial data sets provide a model to describe 
head on black hole collisions. Both initial data sets are 
time symmetric, conformally flat, and the crucial differ- 
ence between each other lies in their topology: 3 asymp- 
totically flat regions connected by 2 throats in the BL 
case, and 2 asymptotically flat regions connected by 2 
throats in the Misner case. For physically equivalent BL 
and Misner data sets it is understood data sets such that: 
(1) the ADM mass of the initial hypersurface measured 
at the reference asymptotic end is the same; and (2) the 
geodesic separation of the innermost minimal surfaces is 
the same. For the case of the Misner data, there is an ex- 
plicit formula for the separation of the minimal surfaces, 
while for the BL data it has to be calculated numerically. 
Cadez [|l4| has given a tabulation of the values of this 
distance for different values of the parameters of the BL 
initial data. Due to the axial symmetry of the Misner 
and BL initial data sets, there is only one non- vanishing 
NP constant for each one of the developments (C?2)- An 
analysis of the NP constants for the BL data has been 
recently given in Eg]. The function W — see formula 
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FIG. 3: The Newman-Penrose constants for physically equiv- 
alent Misner and Brill-Lindquist initial data sets. 



( pli| ) — for the BL data has been calculated there. The 
one for the Misner data can be calculated in a similar 
fashion. A plot of the values of the constant for various 
physically equivalent initial data sets is shown in figure 
Q. It can be seen that for black holes initially very sepa- 
rated the difference between the NP constants of the two 
topologically different initial data set is negligible, and 
it increases as one considers data sets where the separa- 
tion is smaller. For black holes separated a distance of 
L w 5.26(m/2) where m is the total ADM mass of the 
spacetime one has \G 2 (bl) \ 2 /\G 2( m isner) \ 2 ~ °- 66 - The 
physical content of this difference, and the possibility of 
interpreting it as some measure of the (spurious) radia- 
tion contained in the initial data sets is certainly a matter 
of debate. 
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